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Abstract 



Theoretical basics and experimental determinations of the coupling parameter of 
the Strong Interaction, Qs, are reviewed. The world average value of Os, expressed 
at the energy scale of the rest mass of the Z'^ boson, is determined from analyses 
which are based on complete NNLO perturbative QCD. The result is as(Af2o) = 
0.1184 ± 0.0031. No significant deviations or systematic biases of subsamples of 
experimental results are found. From the observed energy dependence of Qs, which is 
CNJ ' in excellent agreement with the expectations of QCD, the number of colour degrees 

' of freedom can be constrained to Nc — 3.03 ± 0.12. 

o 
o 
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1 Preface 



The coupling strength as is the basic free parameter of Quantum Chromodynamics 
(QCD), the theory of the Strong Interaction |Q] which is one of the four fundamental 
(~| , forces of nature. QCD describes the interaction of quarks through the exchange of an 

octet of massless vector gauge bosons, the gluons, using similar concepts as known from 
Quantum Electrodynamics, QED. QCD, however, is more complex than QED because 
^ ' quarks and gluons, the analogues to electrons and photons in QED, are not observed as 

free particles but are confined inside hadrons. 

Confinement implies that the coupling strength as, the analogue to the fine structure 
constant a in QED, becomes large in the regime of large-distance or low- momentum 
transfer interactions|]. Conversely, quarks and gluons are probed to behave like free 
particles, for short time intervalsQ in high-energy or short-distance reactions; they are 
said to be "asymptotically free", i.e. as — > for momentum transfers Q — » c». 

Within QCD, the phenomenology of confinement and of asymptotic freedom is re- 
alized by introducing a new quantum number, called "colour charge" . Quarks carry 
one out of three different colour charges, while hadrons are colourless bound states of 3 
quarks or 3 antiquarks ("baryons"), or of a quark and an anti-quark ("mesons"). Gluons, 
in contrast to photons which do not carry (electrical) charge by themselves, have two 

^ In the world of quantum physics, "large" distances As correspond to As > 1 fm, "low" momentum 
transfers to Q < 1 GeV/c... 

^ ... and "short" time intervals correspond to At < 10~^* s. 
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colour charges. This concept leads to the process of gluon self-interaction, which in turn, 
through the effect of gluon vacuum polarization, gives rise to asymptotic freedom, i.e. 
the decrease of as with increasing momentum transfer. 

As in the case of QED, QCD predicts the energy dependence of ckg, while the actual 
value of as, at a given energy or four momentum transfer scale^ Q, is not predicted but 
must be determined from experiment. 

Determining as at a specific energy scale Q is therefore a fundamental measurement, 
to be compared with measurements of the electromagnetic coupling a, of the elementary 
electric charge, or of the gravitational constant. Testing QCD as such, however, requires 
the measurement of as at least at two different energy scales, and/or at different processes: 
one measurement fixes the free parameter and thus provides accurate predictions for the 
value of as at other energy scales and/or at other processes. 

In general, as can be determined in dynamic particle reactions involving in- or outgo- 
ing quarks and gluons, which manifest themselves as hadrons. Examples of Feynman di- 
agrams describing hadronic final states in deep inelastic lepton-nucleon scattering (DIS), 
electron-positron annihilation (e''"e^), hadron collisions and quarkonia decays are shown 
in figure |^. 




Figure 1: Examples of Feynman diagrams describing hadronic final states in processes 
which are used to measure as. 

In this report, at the turn of the millennium, the current status of measurements of as 
is reviewed. Theoretical basics of ag and QCD are given in Section ||. Measurements of 
as from deep inelastic scattering, from e^e^ annihilation processes, from hadron colliders 
and from heavy quarkonia decays are discussed in Sections 3 to 6, respectively. A global 
summary of these results, a determination of the world average value of as(M2o) and 
quantitative studies of the energy dependence of as are presented in Section |^. Section ^ 
concludes and gives an outlook to future developments. 



2 QCD and as', basic theoretical predictions 

The concepts of QCD are described in many text books and review articles, see e.g. 
references |§, ||] • A brief review of the basics of perturbative QCD and of the coupling 

^Here and in the following, a system of units is utilized where the speed of light and Planck's constant 
are put to unity, c = ?i = 1, such that energies, momenta and masses are all given in units of GeV. 
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strength ag, like the concepts of renormahzation, asymptotic freedom and confinement, 
the A parameter, the treatment of quark masses and thresholds, perturbative predictions 
of physical observables and renormahzation scale dependence, and of nonperturbative 
methods like lattice calculations will be presented in the following subsections. 

2.1 Renormalization 

In quantum field theories like QCD and QED, dimensionless physical quantities TZ can 
be expressed by a perturbation series in powers of the coupling parameter as or a, 
respectively. Consider TZ depending on as and on a single energy scale Q. This scale 
shall be larger than any other relevant, dimensionful parameter such as quark masses. 
In the following, these masses are therefore set to zero. 

When calculating 7^ as a perturbation series in as, ultraviolet divergencies occur. 
Because TZ must retain physical values, these divergencies are removed by a procedure 
called "renormalization" . This introduces a second mass or energy scale, fi, which repre- 
sents the point at which the subtraction to remove the ultraviolet divergencies is actually 
performed. As a consequence of this procedure, TZ and as become functions of the renor- 
malization scale fi. Since TZ is dimensionless, we assume that it only depends on the ratio 
Q"^ I \j? and on the renormalized coupling as(/i^): 



Because the choice of /i is arbitrary, however, TZ cannot depend on /z, for a fixed value 
of the coupling, such that 



where the convention of multiplying the whole equation with [i? is applied in order to 
keep the expression dimensionless. Equation ^ implies that any explicit dependence of 
TZ on ij, must be cancelled by an appropriate /i-dependence of as- It would therefore 
be natural to identify the renormalization scale with the physical energy scale of the 
process, fi^ = Q^, eliminating the uncomfortable presence of a second and unspecified 
scale. In this case, as transforms to the "running coupling constant" cks(Q^), and the 
energy dependence of TZ enters only through the energy dependence of as{Q^). 

2.2 tts and its energy dependence 

While QCD docs not predict the actual size of as at a particular energy scale, its energy 
dependence is precisely determined. If the renormalized coupling as{n'^) can be fixed 
(i.e. measured) at a given scale /i^, QCD definitely predicts the size of as at any other 
energy scale through the renormalization group equation 



TZ = 7^((3^/^^, as); as = as(p.^). 




S/i^ das J 



(1) 




(2) 
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The perturbative expansion of the (3 function, including higher order loop corrections to 
the bare vertices of the theory, is calculated to complete 4-loop approximation M: 



(3) 
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(4) 



and Nf is the number of active quark flavours at the energy scale Q. The numerical 
constants in equation ^ are functions of the group constants Ca and Cf, which for QCD 
— exhibiting SU{i) symmetry — have values of Ca = 3 and Cp = 4/3 (see e.g. [D for 
more details) . (3q and /3i are independent of the renormalization scheme, while all higher 
order /3 coefficients are scheme dependent. 



2.3 Asymptotic freedom and confinement 

A solution of equation ^ in 1-loop approximation, i.e. neglecting [ii and higher order 
terms, is 



2^ 



1 + as(Ai^)/3oln^ 

Apart from giving a relation between as(Q^) and Q!s(/i2), equation |^ also demonstrates 
the property of asymptotic freedom: if becomes large and /3o is positive, i.e. if 
Nf < 17, as((5^) will decrease to zero. 

Likewise, equation]^ indicates that as{Q'^) grows to large values and actually diverges 
to infinity at small Q^: for instance, with as(/^^ = 7Vf|o) — 0.12 and for typical values 
of TV/ = 2 ... 5, as(Q^) exceeds unity for < 0(100 MeV ... 1 GeV). Clearly, this 
is the region where perturbative expansions in as are not meaningful anymore, and we 
may regard energy scales of and Q"^ below the order of 1 GeV as the nonperturbative 
region where confinement sets in, and where equations ^ and ^ cannot be applied. 

Including /3i and higher order terms, similar but more complicated relations for 
as(Q^), as a function of as{n^) and of In % as in equation ^, emerge. They can be 
solved numerically, such that for a given value of as{fi'^), choosing a suitable reference 
scale like the mass of the Z° boson, ^ — Mzo, ctsiQ'^) can be accurately determined at 
any energy scale > 1 GeV^. 
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2.4 The A parameter 



Alternatively, as another parametrization of as{Q^), a parameter called A is introduced 
as a constant of integration of the /3-function, such that, in leading order, 



Equations ^ and ^ are equivalent with each other if 



l/(/3oa.(/x2)) 



Hence, as(/^^) is replaced by a suitable choice of the A parameter, which technically 
is identical to the energy scale Q where as{Q^) diverges to infinity, as{Q^) — > oo for 
q2 _^ ^2 gj^^, ^ numerical example, A « 0.1 GcV for as(A/zo = 91.2 GcV) = 0.12 
and Nf = 5. 

The parametrization of the running coupling as{Q^) with A instead of as(/^^) has 
become a common standard, see e.g. reference Q, and will also be adapted in this 
review. While being a convenient and well-used choice, however, this parametrization 
has several shortcomings which the user should be aware of. 

First, requiring that as{Q'^) must be continuous when crossing a quark thresholc^, 
A actually depends on the number of active quark flavours. Secondly, A depends on 
the renormalization scheme, see e.g. reference In this review, the so-called "modified 
minimal subtraction scheme" (MS) will be adopted, which also has become a common 
standard A will therefore be labelled A^^'* to indicate these peculiarities. 

In complete 4-loop approximation and using the A-parametrization, the running cou- 
pling is thus given by 



1 1 



/3o 



' ^^^^ln3L + ^n^L + 21nL4V3^1nL. 



V ' 2 ^ 2; /32 ^ 2/3o 



(7) 



where L — Q'^/ A^L—.. The first line of equation ^ includes the 1- and the 2-loop coefficients, 
the second line is the 3-loop and the third line is the 4-loop correction, respectively. 

The functional form of as{Q), for the 1-, the 2- and the 3-loop approximation of 
equation 0, each with A = 0.220 GeV, is shown in figure ||(a). As can be seen, there is 
an almost 15% decrease of as when changing from 1-loop to 2-loop approximation, for 
the same value of A. The difference between the 2-loop and the 3-loop prediction is only 
about 1-2%, and is less than 0.01% between the 3-loop and the 4-loop presentation which 
cannot be resolved in the figure. 

^Strictly speaking, physical observables TZ rather than as must be continuous, which may lead to 
small discontinuities in as{Q'^) at quark thresholds in finite order perturbation theory; see section 2.5. 
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C^(4-(oop) _^(n-!oop)\ 

The fractional difference in the energy dependence of ag, — 14-ioav) j for 

n — 1, 2 and 3, is presented in figure ||(b). Here, in contrast to figure ||(a), the values 
of were chosen such that q;s(Mzo) = 0.ff9 in each order, i.e., = 93 MeV (1- 
loop), Aj^ = 239 MeV (2-loop), and Aj^ = 220 MeV (3- and 4-loop). Only the 1-loop 
approximation shows sizeable differences of up to several per cent, in the energy and 
parameter range chosen, while the 2- and 3-loop approximation already reproduce the 
energy dependence of the best, i.e. 4-loop, prediction quite accurately. 



2.5 Quark masses and thresholds 

So far in this discussion, finite quark masses niq were neglected, assuming that both the 
physical and the renormalization scales and /i^ , respectively, are larger than any other 
relevant energy or mass scale involved in the problem. This is, however, not entirely 
correct, since there are several QCD studies and as determinations at energy scales 
around the charm- and bottom-quark masses of about 1.5 and 4.7 GeV, respectively. 

Finite quark masses may have two major effects on actual QCD studies: Firstly, 
quark masses will alter the perturbative predictions of observables TZ. While phase space 
effects which are introduced by massive quarks can often be studied using hadronization 
models and Monte Carlo simulation techniques, explicit quark mass corrections in higher 
than leading perturbative order are available only for jet production and for total 



Secondly, any quark- mass dependence of TZ will add another term A'^f^^^ to 
equation |l|, which leads to energy-dependent, running quark masses, m„{Q^), in a similar 
way as the running coupling as{Q^) was obtained, see e.g. reference 

In addition to these effects, as indirectly also depends on the quark masses, through 
the dependence of the /3 coefficients on the effective number of quarks flavours, TV/, 
with ruq <^ 11. Constructing an effective theory for, say, (iV/-l) quark flavours which 
must be consistent with the Nf quark flavours theory at the heavy quark threshold 
^(^/) ^ 0{mq), results in matching conditions for the as values of the (Nf-l)- and the 
A^/-quark flavours theories . 

In leading and in next-to-leading order, the matching condition is ai^^ — . 
In higher orders and the MS scheme, however, nontrivial matching conditions apply 
Ip^ , ^ 1^. Formally these are, if the energy evolution of is performed in n*'* order, of 
order (n — 1). 

The matching scale fx^-^f ^ can be chosen in terms of the (running) Af S mass mq{iiq), or 
of the constant, so-called pole mass Mq. For both cases, the relevant matching conditions 
in NNLO are given in These expressions have a particluarly simple form for the 
choic^ /i'-^-f^ — mq{mq) or /i'^^J'^ — Mq. In this report, the latter choice will be used 
to perform 3-loop matching at the heavy quark pole masses, in which case the matching 
condition reads, with a = ai^^'/^ ^-i^^ a' = ai^^ ^V""- 

^The results of reference jsj are also valid for other relations between ^^^f^ and rriq or Mq, as e.g. 
^(Nf) _ 2Mq. For 3-loop matching, however, practical difEercnces due to the freedom of this choice are 
negligible. 
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- = 1 + C2 a' + Ca , (8) 
a 

where C2 = -0.291667 and C3 = -5.32389 + {Nf - 1) • 0.26247 g. 

The 4- loop prediction for the running as, using equation ^ with A^^"^-* = 220 MeV 

and 3-loop matching at the charm- and bottom-quark pole masses, /xi^^"*'' — — 
1.5 GeV and = A'h ^ 4.7 GeV, is illustrated in figure |a (full hue). SmaU 

discontinuities at the quark thresholds can be seen, such that ai^' < ai'^^'^ by about 
2 per mille at the bottom- and about 1 per cent at the charm-quark threshold. The 
corresponding values of Aj^ are A^^"''^ = 305 MeV and = 346 MeV p4 |. 

Comparison with the 4-loop prediction, without applying threshold matching and for 
A.^"^"* — 220 MeV and Nf — 5 throughout (dashed line) demonstrates that, in spite of 
the discontinuities, the matched calculation shows a steeper rise towards smaller energies 
because of the larger values of A^"**^ and a!^~^\ 

MS MS 

The size of discontinuities and the changes of slopes are more clearly demonstrated 
in figure where the fractional difference between the two curves from figure |^a, i.e. 
between the matched and the unmatched calculation, is presented. Note that the step 
function of as is not an effect which can be measured — the steps are artifacts of the 
truncated perturbation theory and the requirement that predictions for observables at 
energy scales around the matching point must be consistent and independent of the two 
possible choices of (neighbouring) values of Nf. 



2.6 Perturbative predictions of physical quantities 

In practice, as is not an 'observable' by itself. Values of as{n^) are determined from mea- 
surements of observables TZ for which perturbative predictions exist. These are usually 
given by a power series in as(/i^), like 

n{Q^)^ Pi Y^Rna: 

n 

= P/(i?o + i?ias(/i') + ^2(QVAi'K'(A'') + -) , (9) 

where i?„ are the nth order coefficients of the perturbation series and PiRo denotes the 
lowest-order value of TZ. 

For processes which involve gluons already in lowest order perturbation theory. Pi 
itself may include (powers of) ag. For instance, this happens in case of the hadronic 
decay width of heavy Quarkonia, r(T — > ggg — > hadrons) for which Pi oc a^. If no 
gluons are involved in lowest order, as e.g. in e"'"e~ ^ 99 — > hadrons or in deep inelastic 
scattering processes, PiRq is a constant and the usual choice is P; = 1. Rq is called the 
lowest order coefficient and Ri is the leading order (LO) coefficient. Following this naming 
convention, R2 is the next-to-leading order (NLO) and R3 is the next-to-next-to-leading 
order (NNLO) coefficient. 

QCD calculations in NLO perturbation theory are available for many observables 
TZ in high energy particle reactions; calculations including the complete NNLO are 



7 



available for some totally inclusive quantities, like the total hadronic cross section in 
e"'"e~ — > hadrons, moments and sum rules of structure functions in deep inelastic scat- 
tering processes and the hadronic decay width of the r lepton. The complicated nature 
of QCD, due to the process of gluon self-coupling and the resulting large number of 
Feynman diagrams in higher orders of perturbation theory, so far limited the number of 
QCD calculations in complete NNLO. 

An alternative approach to calculating higher order corrections is based on the re- 
summation of leading logarithms which arise from soft and collinear singularities in gluon 
emission p^ . In such a case, the effective expansion parameter is ofgL^ rather than ag, 
where L = ln(l/7?.) and TZ is some generic observable which tends to zero in lowest order. 
For small values of TZ, asL/^ becomes large, and therefore these terms should be known 
to all orders in as if a reliable prediction of TZ is to be obtained. For certain observables 
it has proved possible to sum up both the leading and next-to-leading logarithms, which 
is referred to as the 'Next-to-Leading Log Approximation' or NLLA. 

For observables TZ for which resummation is possible, the cumulative cross-section 
E(7?.) may be written as 

E(7^)=/ -^c^7^ = C(as)exp[G(a„L)]-|-D(as,7^) , (10) 
Jo f "'^ 

where D(as,TZ) is a remainder function which should vanish as 7?. ^ 0, and 

oo 

C{as) = 1 + ^ Cna: (11) 

oo n+1 

G(a„L) = Y.Y1 Gnma^L^ 

n—1 m— 1 

= Lgi{asL) + g2{asL) + asgsiasL) + algi{asL) • ■ • , 

with ots = (as/27r). The functions Lgi and g2 represent the sums of the leading and 
next-to-leading logarithms respectively, to all orders in a^- 

Resummation of leading and next-to-leading logarithms provides predictions which 
include leading terms of all perturbative orders; however, they are not complete in the 
sense of fixed-order predictions because the latter also include sub-leading logarithms 
and non-logarithmic terms. This is illustrated in Table |l], where the NLLA calculations 
provide the sum of all terms of the first two columns, while NLO {0{a1)) calculations 
yield the sums of the terms in the first two rows. 

If the maximum available information from both resummed and fixed order calcula- 
tions shall be utilized, they must be added such that double counting of terms which are 
common to both these calculations is avoided. This can be achieved involving several 
(approximate) 'matching schemes'. For instance, terms to 0{a1) in the NLLA expres- 
sion for lnS(7?.) (c.f. table ^ can be removed and replaced by the full expression in 
exact 0{al). This procedure is called 'Ini? matching'; it results in 'resummed ©(a^)' 
predictions which should be superior to the fixed 0{a1) case because more higher order 
correction terms are included. 

The perturbative order up to which QCD predictions for different processes and 
observables are available, is indicated in Table o, at the end of this review. 
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Table 1: Decomposition of the cumulative cross-section, lnS(7^), in powers of = 
(as/27r) and L = ln(l/7e). 
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2.7 Renormalization scale dependence - infinite discomfort in fi- 
nite order. 

The principal independence of a physical observable TZ from the choice of the renormal- 
ization scale /i was expressed in equation Replacing by Q;s(/i^), using equation ^, 
and inserting the perturbative expansion of TZ (equation ^ into equation |l| results in 

i?i/3o 

-f 2i?2/3o] 

(12) 

Solving this relation requires that the coefficients of a"(/i^) vanish for each order n. With 
an appropriate choice of integration limits one thus obtains 

i?o = const. , 
i?i = const. , 

i?2 =i?2(l)-/3oi?iln^ , 

^3 = ^3(1) - [2i?2(l)/3o + RiPi] In ^ + RiPl In' ^ (13) 

as a solution of equation |lj. 

In other words, invariance of the complete perturbation series against the choice of the 
renormalization scale /i^ implies that the coefficients except Rq and explicitly 
depend on /i^. In infinite order, the renormalization scale dependence of as and of 
the coefficients i?„ cancel; in any finite (truncated) order, however, the cancellation 
is not perfect, such that all realistic perturbative QCD predictions include an explicit 
dependence on the choice of the renormalization scale. 



= ^l■ 



,dRo 
diJ? 



, dRi 



Oiat) 



, dR2 
i dRs 
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This renormalization scale dependence is most pronounced in leading order QCD 
because Ri does not depend on /i and thus, there is no cancellation of the (logarithmic) 
scale dependence of Q!s(/^^) at all. Only in next-to- leading and higher orders, the scale 
dependence of the coefficients i?„, for n > 2, partly cancels that of as(/i^). In general, the 
degree of cancellation improves with the inclusion of higher orders in the perturbation 
series of 7?.. 

A practical example of the scale dependence of as, determined from the measured 
value of the scaled hadronic decay width of the Z" boson jl^ , 

^z-^ir^^^^^ 20.768 ±0.0024, (14) 
1 (Z'J leptons) 

is shown in figure ^. QCD predictions of Rz are available in complete NNLO p^ . 
Including further corrections like quark mass and non-factorizable electroweak and QCD 
effects, these predictions can be parametrised to 



Rz = 19.934 



1.045^^ + 0.94 





2 

- 15 




3' 


TT 




TT 





(15) 



for ^ = Mza; see section 4.3 for more details. In order to demonstrate the scale depen- 
dence of as, as(/i) is determined from equations [l4| and [T5| in LO, in NLO and in NNLO, 
as a function of the scaled renormalization scale, ~ /i/AJ^o, using the appropriate, 
scale dependent QCD coefficients i?„ according to equation |l^. as(A/zo) is then calcu- 
lated from as(/i) using equation in 1-loop expansion for the LO, in 2-loop for the NLO 
and in 3-loop for the NNLO case. 

The resulting scale dependence of as {Mzo ) is displayed in figure ^(a) ; the scale depen- 
dence of the QCD coefficients Ri , i?2 and R3 is given in figure ^(b) . The scale dependence 
of as is visibly reduced if higher orders are included: limiting renormalization scales to 
factors from 1/5 to 5 of the "physical" scale Q = Mzo , results in changes of as {Mzo ) of 
about 40% in LO, 8% in NLO and 2.5% in NNLO. Beyond these hmits, however, the 
scale dependence even in NNLO QCD diverges, such that it seems not to be meaningful 
to consider yet larger ranges of scale factors - although there is no rule, from first 
principles, which limits scale variations to a definite range. 

There are several proposals to optimize or to fix the renormalization scale in finite 
order perturbation theory. An intuitive approach is to apply the general scale insensi- 
tivity requirement, equation |l|, which is strictly valid in infinite order, to calculations in 
any finite (beyond leading) order. This is Stevenson's principle of minimal sensitivity 
(PMS)||l^, where, for each observable TZ, the optimal scale is defined by dTZ/d^^ — 0. 
The PMS solution for Rz in NLO is marked in figure ||(a) [filled circle] . 

In the "effective charge" method (EC) of Grunberg jl9j, all higher (i.e. beyond lead- 
ing) order coefficients Rn vanish. The EC solution in NLO is also marked in figure ^(a) 
[open circle]. 

Brodsky, Lepage and Mackenzie (BLM) |2^ propose to choose /i such that the NLO 
coefficient R2 is independent of Nf, the number of active quark flavours. The BLM 
solution is marked in figure ^(a) [open rhombus] . 
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PMS, EC and BLM all provide unambiguous solutions in NLO; in the case of Rz, 
they are very close to each other but the NNLO prediction does not exactly match these 
results: the value of as from "optimized" NLO can only be reached in NNLO at very 
small values of a;^;. 

While in NLO QCD it is sufficient to vary the renormalization scale in order to assess 
the full spectrum of variation offered by the renormalization procedure, in NNLO both 
scale- and scheme-dependences must be taken into account. PMS, EC and BLM can 
then be generalized to NNLO and beyond, however do not lead to unambiguous results 
anymore. The whole issue of theoretical scale optimizations was and still is vividly 
discussed in the literature, see e.g. plf . 

Another approach is motivated by an experimental point of view: ideally, the optimal 
choice of /X in a given order of perturbation theory should reproduce or, at least, closely 
approach the unknown all-order result. Measurements of observables should inherently 
include all orders of perturbation theory; therefore it may be possible to extract the 
optimal choice of /i in a given, finite order calculation from fits to the experimental data. 
Indeed, experimental scale optimization is possible in cases where as is determined from 
differential distributions of observables, and often leads to a remarkable consistency 
between results obtained from different observables |2^, |2^, [2^, . 

So far, none of the methods described above was generally accepted as the preferred 
method to optimize or fix renormalistaion scales in finite order calculations. Instead, 
in experimental determinations of ag, systematic uncertainties due to unknown higher 
order contributions are usually defined by allowing the renormalization scale fi to vary 
in "reasonable" ranges. 

Unfortunately, there is no common agreement on what should be a reasonable range, 
and the significance and interpretation of such procedures is unclear. Furthermore, as 
indicated above, scale changes alone are not sufficient to investigate higher order uncer- 
tainties in NNLO predictions, where renormalization scheme dependences should also be 
accounted for. 

2.8 Lattice QCD 

At large distances or low momentum transfers, as becomes large and application of 
perturbation theory becomes inappropriate. In this regime, nonperturbative methods 
must be used to describe processes of the Strong Interaction. Lattice QCD is one of 
the most developed nonperturbative methods which is used to calculate, for instance, 
hadron masses, hadron mass splittings and QCD matrix elements. In Lattice QCD, field 
operators are applied on a discrete, 4-dimensional Euclidean space-time of hypercubes 
with side length a. 

Energy levels of heavy quarkonia systems (QQ) calculated using lattice QCD depend 
on the heavy quark mass Mq and on as. Comparison with measurements of quarkonia 
mass splittings then allows to determine ag, which at this point is not given in the MS 
renormalization scheme but is based on other suitable definitions of as on the lattice. 
Conversion from the lattice to the MS coupling at high energies can be done using an 
expansion in second or in third order perturbation theory. A review of methods to 
determine as from lattice investigations is given in [p6| 
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So far, calculations exist which either neglect light quark loop contributions (Nf—0; 
"quenched approximation") or which include two light quark flavours (Nf=2); the latter 
allow extrapolation to Nf=S. Uncertainties in this extrapolation, the limited order to 
which the conversion to the MS scheme is known, limited Monte-Carlo statistics and 
corrections for light quark masses lead to theoretical imcertainties of final (MS) values 
of as(Mzo). 



3 Results from deep inelasting scattering processes 

Measurements of the violation of Bjorken scaling in deep inelastic lepton-nucleon scatter- 
ing belong to the earliest methods to determine as- The first significant determinations 
of cks, i.e. those based at least on next-to-leading order (NLO) perturbative QCD pre- 
diction, started to emerge in 1979 p7| ]. 

Today, a large number of results is available, from data in the energy (Q^) range of 
a few to several thousand GeV^, using electron-, muon- and neutrino-beams on various 
fixed target materials (H, D, C, Fe and CaCOa), as well as electron-proton or positron- 
proton colliding beams at HERA. In addition to scaling violations of structure functions, 
as is also determined from moments of structure functions, from QCD sum rules and — 
as in e^e~ annihilation — from hadronic jet production and event shapes. 



3.1 Scaling violations and moments of structure functions 

Cross sections of physical processes in lepton-nucleon scattering and in hadron-hadron 
collisions depend on the quark- and gluon-densities in the nucleon. Assuming factor- 
ization between short-distance, hard scattering processes which can be calculated using 
QCD perturbation theory, and low-energy or long-range processes which are not accessi- 
ble by perturbative methods, such cross sections are parametrized by a set of structure 
functions Fi [i— 1,2,3). The transition between the long- and the short-range regimes is 
defined by an arbitrary factorization scale /i / , which — in general — is independent from 
the renormalization scale but has similar features as the latter: the higher order coef- 
ficients of the perturbative QCD series for physical cross sections depend on fif in such 
a way that the cross section to all orders must be independent of /i/, i.e. da/dfj,f = 0. 
To simplify application of theory to experimental measurements, the assumption fif — fi 
is usually made, with /i = Q as the standard choice of scales. 

In the naive quark-parton model, i.e. neglecting gluons and QCD, the differential 
cross sections for electromagnetic charged lepton (electron or muon) or weak neutrino- 
proton scattering off an unpolarized proton target are written 

j2 em 1 + n — 1/]^ M 

^ - ^^2.F, + (1 - y)iFr ~ 2.Fr) ~ ^.yF^ ; 

|£ ^ (i - . - ^^y) n + y^^n + . (i - i.) , m 

where x — 2M(e-E') momentum fraction of the nucleon carried by the struck 

parton, y = 1 — E' / E, is the negative quadratic momentum transfer in the scattering 
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process, and M, E and E' are the mass of the proton and the lepton energies before and 
after the scattering, respectively, in the rest frame of the proton. In the quark-parton 
model, these structure functions consist of combinations of the quark- and antiquark 
densities q{x) and q{x) for both valence- (u,d) and sea-quarks (s,c). 

In QCD, the gluon content of the proton as well as higher order diagrams describ- 
ing photon/gluon scattering, qq, and gluon radiation off quarks must be taken 
into account. Quark- and gluon-densities, the structure functions Fi and physical cross 
sections become energy (Q^) dependent. QCD thus predicts, departing from the naive 
quark-parton model, scaling violations in physical cross sections, which are associated 
with the radiation of gluons. While perturbative QCD cannot predict the functional 
form of parton densities and structure functions, their energy evolution is described by 
the so-called DGLAP equations fH, H]. 

Structure functions contain, apart from terms whose energy dependence is given by 
perturbative QCD, which so far is known up to complete NLO js^, so-called "higher 
twist" contributions (HT). The leading higher twist terms are proportional to l/Q^; 
they are numerically important at low < ©(few GeV^) and at very large a; ~ 1. 

Historically, precise results of as, from the logarithmic slopes of F2{x,Q^) and in 
NLO QCD, were obtained from a combined analysis ||3^ of SLAC and BCDMS data, in 
a range from 0.5 to 260 GeV^, giving Q!s(Mzo) = 0.113 ± 0.005. The error includes 
theoretical (scale) uncertainties of ±0.004 |^^. Higher twist terms as well as the gluon 
distribution were simultaneously determined in this analysis. 

The CCFR colaboration obtained as from i^-nucleon scattering and a fit to the non- 
singlet structure function F^{x, Q^) | |3^ , which is independent of the poorly known gluon 
distribution. In order to increase statistical precision, F^ was substituted by F2 at 
large x where gluons do not contribute much, giving as(A/zo) = 0.111 ± 0.002(stat.) ± 
0.003(sys.) ±0.004(theo.). 

These two results were, for several years, the most significant from deep inelastic 
lepton-nucleon scattering. Because they were numerically smaller than typical values 
obtained from e+e" annihilation, as(Afzo) ^ 0.120 (see Section 4), speculations about 
possible reasons for such differences arose. These speculations came to a halt, at least 
partly, when the CCFR collaboration corrected their previous result — due to a new 
energy calibration of the detector — to as(Mzo) = 0.119 ± 0.005 [Q. Further results 
from F2, in a study |^ of the HERA data at small x and < 100 GeV^, based on 
NLO QCD including summations of all leading and subleading logarithms of and 
1/x, lead to as(Mzo) = 0.120 ± 0.005 ± 0.009, where the first error is experimental 
and the second theoretical. Finally, a reanalysis of the SLAC, BCDMS and NMC data 
on F2, taking proper account of point-to-point correlations, resulted in an increase of 
as(Mzo) from 0.113 to 0.118 |Q. A continuation of this analysis, including recent HERA 
data and careful studies of the effects of higher twist terms, was reported to result in 
a,(Mzo) = 0.116 ±0.003 [|7|. 

Significant progress in this field was achieved with the availability of NNLO QCD 
predictions for the non-singlet structure function F3 as well as for moments of i^2- A new 
analysis of the CCFR data, now in NNLO QCD and from 2:^3 alone, resulted in 

as(A/zo) = 0.118 ± 0.002(stat.) ± 0.005(sys.) ± 0.003(thco.) ; 
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see figure |g. This value is taken as the currently most significant result from struc- 
ture functions in neutrino-nucleon deep inelastic scattering, and is included in the final 
summary, see Table ^. Another recent study was based on Bernstein polynomials and 
moments of F2 from electron- and muon-scattering data, including fixed target as well 
as HERA coUiding beam data in the range of 2.5 to 230 GeV^. The result |Q was 

a^iMzo) = 0.1172 ± 0.0017 ± 0.0017 , 

where the first error is experimental and the second theoretical. The theoretical error 
includes higher twist effects and an estimate of the NNNLO corrections. This result is 
taken as the final value from F2 in DIS and is added to the summary in section 7. 

From scaling violations of polarized structure functions, based on data from SLAG 
and SMC, Ofg was determined in NLO QCD, resulting in 

a,(A/zo) = 0.120iHo5(exp.)-o:oo6(theo.) , 
which is also added to the final summary. 



3.2 Sum Rules 

Perturbative corrections to two inclusive measurements, namely the Gross-Llewellyn 
Smith (GLS) sum rule for deep inelastic neutrino scattering and the Bjorken sum 
rule ||42| for polarized structure functions, have been determined to complete NNLO 
QGD |3[ ||. The GLS sum rule, 

GLS= / Fs{x,Q^)dx = 3(1- — + ...), (17) 
Jo 

when fitted to data of the GGFR collaboration |^ at = 3 GeV^, resulted in p3|] 
as(l-73 GeV) — 0.32 ± 0.05. A recent update of the data and analysis from CCFR ||46|] , 
in the range from 1 to 15.5 GeV^, gave 

as(L73 GeV) = 0.28 ± 0.035(stat.) ± 0.050(sys.)+°;°3o(theo.) . 

The systematic uncertainty is dominated by the extrapolation of the GLS integral to the 
regions x < 0.01 where no measurements exist, and to a; > 0.5 which is substituted by 
F2 from SLAG data, see figure ^. The theoretical error is dominated by uncertainties in 
the higher twist corrections. 

The Bjorken polarized sum rule determines that 

J\glix)~g-ix)]dx^^ 

where the polarized structure functions 51(2;) for protons (superscript p) and neutrons 
(n) are derived from the difference of cross sections for parallel and antiparallel polarized 
targets and (/x or e) beamsjOA and gv are the constants of the neutron weak decay, 
g^/gy = -1.26. A study of the GERN SMG [|| and the SLAG E142 || data 
obtained 

as(1.58 GeV) = 0.375io;°8i. 

No explicit corrections for nonperturbative higher twist where applied to derive this 
result, however an estimate of the size of the 0{ag) terms pO] was taken into account. 



gA 

gv 



(18) 
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3.3 Jet Production and Event Shapes 



Observables parametrizing hadronic event shapes and jet production rates are the clas- 
sical inputs for ctg studies in e''"e~ annihilation. In recent years, these observables were 
also studied in high energetic electron- and positron-proton coUisons at HERA, leading 
to significant determinations of ag- 
in detail, inclusive as well as difi^erential jet production rates were studied in the 
energy range of ~ 10 up to 10000 GeV^ |Q, using variants of jet algorithms from 
e"'"e~ annihilation which are decribed, in more detail, in section 4. In leading order as, 
2-1-1 jet events in deep inelastic ep scattering arise from photon-gluon fusion and from 
QCD compton processes. The term '2-1-1 jet' denotes events where two resolved jets 
can be identified, in addition to the beam jet from the remnants of the incoming proton. 
Previous NLO QCD predictions which were shown to be imprecise are now replaced 
by more recent calculations . 

The results of from jet production at HERA |^ can be summarized to 

as(Mzo) = 0.118 ±0.002(stat.)±0.008(sys.)±0.007(theo.) . 

The systematic error contains uncertainties from using different jet algorithms and hadro- 
nization models, and the theoretical error is dominated by uncertainties from structure 
functions and from scale variations. 

So far, studies of as from hadronic event shape distributions at HERA are based 
on fits to the energy (Q^) evolution of mean values of shape distributions, using NLO 
QCD predictions together with parametrizations of power corrections l/Q^ These 
methods do not yet lead to consistent results of as from different shape observables Q 
— at least not to the degree of precision which is obtained from "classical" studies of as 
from event shape distributions, as e.g. in e"'"e~ annihilation; see Section 4.1 for further 
discussion. Awaiting further progress in theoretical understanding of power corrections 
to event shapes, the as results of these studies will not yet be included in the final 
summary of as . 



4 Results from e+e annihilation 

In e"'"e~ annihilation reactions, as is classically determined from hadronic event shapes, 
jet production rates and energy correlations for which complete NLO QCD calculations 
are available ^ |5^. The first data analysis of this type, from jet rates observed 
at the PETRA e+e~ collider, emerged in 1982 Reviews of early results on as from 
e+e^ annihilations can be found e.g. in references |6l|| . 

More recently, this field was extensively covered by the experiments at the LEP and 
SLC colliders, see e.g. |6|, |^ ^ for reviews of the first few years of LEP op- 
eration. Using data from LEP and from earlier e"'"e~ colliders at lower cm. energies, 
as can be determined from scaling violations of fragmentation functions. Most impor- 
tantly, precise determinations of as are now obtained from the hadronic partial width of 

decays, from overall fits to precision electroweak data and from hadronic branching 
fractions of t lepton decays, which are calculated in complete NNLO QCD. 
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4.1 Event shapes, jet rates and energy correlations 



4.1.1 Observables in (resummed) NLO QCD 

Hadronic event shape variables, jet rates and energy correlations are tools to study both 
the amount of gluon radiation and details of the hadronization process. The definitions of 
observables which are applied to hadronic final states of e"'"e~ annihilations, like Thrust, 
Thrust Major and Minor, Oblateness, jet masses, the jet broadening measures, energy 
correlations and jet production rates, are summarized elsewhere; see e.g. ^ |69|| . 
Only two of these, the Thrust observable and the JADE jet algorithm shall be described 
here in some more detail: 

The Thrust T jTOt is the normalized sum of the momentum components piH of all 
particles z of a given event along a specific axis n; the axis is chosen such that T is 
maximized: 



Thrust assumes values of T = 1 for perfectly aligned momentum vectors, i.e. for an ideal 
event e"''e~ — > qq with two narrow back-to-back jets, down to T = 0.5 for a completely 
spherical distribution of momentum vectors in the limit of events with many gluons 
radiated off the initial qq pair. 

Within the JADE jet algorithm |Q, the scaled pair mass of two resolvable jets i and 
j, Uij = Mfj/E^^^, is required to exceed a threshold value ycut, where Eyis is the sum 
of the measured energies of all particles of an event. In a recursive process, the pair of 
particles or clusters of particles n and m with the smallest value of Unm is replaced by 
(or "recombined" into) a single jet or cluster k with four-momentum pk = Pn +Pm, as 
long as ynm < Vent- The procedure is repeated until all pair masses yij are larger than 
the jet resolution parameter ?/cut, and the remaining clusters of particles are called jets. 

Several jet recombination schemes and definitions of Mij exist |7|] ; the original 

JADE scheme with = 2 ■ ■ Ej ■ {1 — cos 9ij), where Ei and Ej are th e energ ies of 
the particles and 9ij is the angle between them, and the "Durham" scheme |6^ with 
M^j = 2 ■ mm{Ef, Ej) ■ (1 — cos 9ij), were most widely used at LEP, due to their superior 
features like small sensitivity to hadronization and particle mass effects | |6^ . 

QCD predictions for hadronic event shapes, of jet production rates and of energy 
correlations are available in complete NLO |5^, ^ . Differential distributions of such 
observables y are parametrized as 

= Riiy)a,ifi^) + R2{y,Qy^I^)aU^^^) , (19) 

CTo dy 

where ctq is the total hadronic cross section in leading order. In addition, for some of 
the observables, resummation of the leading and next-to-leading logarithms (NLLA) is 
available which can be matched with the NLO expressions (resummed NLO); see 
section 2.6. 

In a typical study, measured distributions are corrected for effects of limited detector 
acceptance and resolution. Fits of QCD predictions to the data are performed after 
applying hadronization corrections, as predicted from hadronization models, to the data 
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or to the theoretical predictions. Systematic uncertainties due to these corrections, from 
renormalization scale variations and from different matching procedures, are studied and 
included in the overall error of the fitted value of as. No common agreement exists, 
however, about the range of variations used to estimate these uncertainties. 

Early determinations of ag from event shapes, jet rates and energy correlations ob- 
tained from experiments at the PETRA and PEP colliders were summarized to as = 
0.14 ± 0.02 at Q2 = ^ (34 GeV)^ |6^, |6|, in NLO QCD. The results which led 
to this average did not include estimates of theoretical uncertainties; however the er- 
ror was determined from the scatter of the single results which were based on different 
observables, and thus gave a first estimate of theoretical uncertainties. 



4.1.2 Results in resummed NLO 

A reanalysis of PETRA data, using refined analysis techniques, resummed NLO QCD 
calculations, modern model calculations and including new observables, quite similar to 
recent analyses at LEP, resulted in 

as (22 GeV) = 0.161 ± 0.009±+°;[5J^ , 
as (35 GeV) = 0.145 ± 0.002+°:°j7 , and 



as (44 GeV) = 0.139 ±0.004 



1.007 ' 
+0.010 
-0.007' 



where the first errors are statistical and experimental systematic, added in quadrature, 
and the second are theoretical uncertainties. These values are taken as the final results 
in the PETRA energy range and are included in the summary table ^. Similar analyses 
from experiments at the TRISTAN collider, around cm. energies of 58 GeV, gave ||7^ 

as (58 GeV) = 0.132 ±0.008 , 

which is also included in the summary. 

At the LEP collider, all four experiments (ALEPH, DELPHI, L3 and OPAL) have 
contributed a multitude of studies based on the high statistics data samples around 
the Z" resonance (LEP-I; Ecm ^ 91.2 GeV) and at the higher energies of the LEP- 
II running phase {Ecm ~ 133, 161, 172, 183 and 189 GeV). The SLD experiment at 
the SLAG Linear Collider (SLG) contributed similar studies at Ecm ~ M^o. As an 
example of the precise description of data by the CQD fits, the Thrust distributions as 
measured by ALEPH at LEP-I and LEP-II energies are shown in figure 0, together with 
the corresponding fits of the resummed NLO QCD calculations. 

The most current results from LEP and SLC, based on resummed NLO QCD cal- 
culations, are taken from |7^, Q |7^, |o) and are combined to obtain one single 
value of tts at each cm. energy. This is done by calculating weighted averages of the 
results quoted by each experiment, with the inverse quadratic total error taken as the 
weight. The experimental uncertainties are combined assuming that no correlation ex- 
ists between the experiments. Theoretical uncertainties are assumed to be common to 
all experiments, however different methods and definitions were used to estimate these 
in each case. Therefore, a linear average of the theoretical uncertainties quoted by the 
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experiments is used to define the theoretical error of the combined resuhs. This gives 



as (91.2 GeV) = 0.121 ± 0.001 ± 0.006 , 
as (133 GeV) = 0.113 ± 0.003 ± 0.006 , 
(161 GeV) = 0.109 ± 0.004 ± 0.005 , 
as (172 GeV) = 0.104 ± 0.004 ± 0.005 , 
as (183 GeV) = 0.109 ± 0.002 ± 0.004 , 
as (189 GeV) = 0.110 ± 0.001 ± 0.004 , 

where the first errors are experimental and the second theoretical. These results are 
included in the final summary table 0. Note, however, that some of the averages at 



LEP-II energies contain results which are still preliminary, see im, [Tl, [78, 79 



The energy dependence of as from studies of event shape observables is clearly seen 
in the results presented above. This is also demonstrated in a dedicated analysis of the 
L3 Collaboration |8^ which includes a study of radiative events recorded at LEP-I. Such 
events are effectively produced at reduced hadronic centre of mass energies; they can 
therefore be used to determine as at energy scales below the nominal collider energy. 
The statistical precision at these reduced energies is, however, limited. 

The results of the L3 study are shown in figure ||, demonstrating the agreement of the 
data with the QCD expectation of a running coupling. Note that, in order to judge the 
significance for the running, only the innermost, experimental errors must be taken into 
account, since theoretical uncertainties are highly correlated between the different data 
points. The running of as was also demonstrated in a recent study of jet production rates 
at PETRA and at LEP energies, using the same consistent analysis method and data 
from two experiments comprising similar detector techniques, JADE and OPAL 

Before the advent of resummed QCD calculations for event shape observables, de- 
terminations of as were done in pure NLO, see e.g. ||2^, Studies of theoretical 
uncertainties typically included scale variations from Xfj_ — 1 down to the best fit values 
of x^, from two-parameter fits, or the optimized scales given by the methods discussed 
in Section 2.7 . These ranges were quite large, down to ~ 0.01, which lead to rather 
large scale uncertainties in case of some of the observables. Resummed NLO results, 
however, had smaller dependencies on scale variations, as theoretically expected, which 
is why they developed as a standard in as determinations from event shape studies. 



4.1.3 Pure NLO results 

Fits of pure NLO QCD to experimental event shape and jet rate distributions are known 
to provide a consistent description of the data if - in addition to as - also the renormal- 
ization scale is treated as a free parameter of the fit j2^ . While the fitted scale factor a;^ 
assumes different values for different observables, indicating that the amount of unknown, 
higher order contributions is different in each case, the resulting value of as appears to 
be "universal" and in much better agreement than applying a single and fixed choice of 
scale H ||. 

These findings are corroborated in a recent re-analysis of DELPHI data at LEP-I ||8^ , 
which also demonstrates that experimentally optimized NLO calculations can provide a 
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more consistent description of the data than resummed NLO at fixed scale — 1. 
Defining a scale range, for each observable, of a factor of 2 around the experimental fit 
value of x^, results in as{Mzo) = 0.1174 ± 0.0026, which includes both experimental and 
theoretical uncertainties. This method of optimizing renormalization scales and defining 
the size of the remaining theoretical uncertainties seems logical from an experimental 
point of view, see the discussion in Section 2.7, but carries intrinsic theoretical problems: 

Firstly, the resulting scale factors are rather small, down to 0(0.01) or smaller, such 
that large logarithms appear in the NLO coefficients, see equation 13. Secondly, there 
is no common agreement about the significance of one overall scale for each distribu- 
tion because theoretical scale optimizations predict to depend on the value of the 
observable itself. The issue of experimental scale optimization therefore requires further 
investigation and discussion. 



4.1.4 Power corrections 



In the past few years, analytical approaches were pursued to approximate nonpertur- 
bative hadronization effects by means of perturbative methods, introducing a universal, 
non-perturbative parameter 



1 f^' 

aoip-i) = — / dfc as{k) 
M/ Jo 



to parametrize the unknown behaviour of as{Q) below a certain infrared matching scale 
III . Divergent soft gluon contributions to the perturbative predictions of event shapes 
are removed and, as a consequence of these techniques, lead to corrections which are pro- 
portional to powers of l/Q. Power corrections are regarded as an alternative approach to 
describe hadronization effects on event shape distributions, instead of using phenomeno- 
logical hadronization models as in the studies discussed in the previous sections. 

The energy dependence of event shape distributions and of their integrated mean 
values, in the cm. energy range from 14 to 183 GeV, were analysed and compared 
with perturbative QCD calculations plus added power corrections |74, 8^, |7^, |7^. Two- 



parameter fits of ao and as{M-zfi) to the data provide a consistent description of the 
shapes and the energy evolution of the data; however the scatter between results from 
different shape observables is larger than the typical uncertainties would suggest. The 
non-perturbative parameter ap turns out to be universal for all studied observables within 
about 20%, and values of as{Mzo) seem to be close to but systematically smaller than 
those from fits using conventional methods of hadronization corrections [p5| . 

Power corrections, as an analytical ansatz to compute nonpcrturbative QCD cor- 
rections to event shape observables, offer a rather promising alternative to the use of 
phenomenological hadronization models in determinations of cks; however more experi- 
ence, confidence and further developments are necessary before actual fit results of as 
can supersede the results from "classical" event shape analyses described above. 



4.2 cts from scaling violations of fragmentation functions 

The total cross section of e+e~ annihilations into charged hadrons /i, e+e^ h + X , 
can be factorized into a perturbative and a nonpcrturbative regime, describing processes 
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of hard gluon radiation and of hadronization, respectively. The perturbative part can 
be calculated in terms of so-called coefficient functions, while the non-perturbative part 
is parametrized by phenomenological fragmentations functions. The separation of these 
two regimes is performed at an energy scale ///, the factorisation scale, which — as in 
DIS — enters as a second, arbitrary energy scale, in addition to the renormalization scale 
jji. According to the factorization theorem, the cross section can then be written js^ : 

^(e+c- ^ + X) = ^ y ' ^ (z, a,{t?), ^i^lQ") A M/) , (20) 

where Ci are the coefficient functions for creating a parton with flavour i and momentum 
fraction z = Pparton/Pbeam, Di represent the probability that parton i fragments into a 
hadron h with momentum fraction x/ and x — Phadron/Pbeam- 

In lowest order, the coefficient functions Ci are given by the electroweak couplings 
for quarks; they vanish for gluons. Higher order QCD corrections apply to the quark 
coefficient functions and lead to finite functions for gluons, too. These corrections are 
known up to complete NLO pTf , however only the LO corrections were available and 
used in currently existing analyses. 

The fragmentation functions Di are not given by perturbation theory; however — as 
in the case of parton densities and structure functions (c.f. Section 3.1) — their energy 
dependence is described by QCD via complicated integro-differential equations and the 
DGLAP evolution functions H, H]. 

The ALEPH Q and the DELPHI ^ collaborations at LEP have both published 
detailed analyses of the charged hadron fragmentation functions, using data samples from 
the PETRA, PEP and LEP colliders spanning cm. energy ranges from 14 to 91.2 GeV. 
At LEP, ALEPH and DELPHI extracted differential x distributions separately for ini- 
tial b-, c- and uds-quark events, and obtained the corresponding gluon-jet distribution 
from tagged gluon jet event samples. Assuming a certain parametrization of the "bare" 
fragmentation functions at a "starting energy scale" Qq, as{Mia), the free parameters 
of the starting fragmentation functions for gluons and the different quark flavours, as 
well as a parameter allowing for additional, nonperturbative power law corrections, are 
determined in flts to all experimental x-distributions. 

Combining both the ALEPH and the DELPHI results, which are based on comparable 
data sets and analysis strategies, leads to 

as{Mzo) = 0.125to;go? (exp.) ±0.009 (theo.) 

where the theoretical uncertainty includes variations of both the factorization and the 
renormalization scalesn. 



4.3 as from total hadronic cross sections 

Because of its inclusive nature, the total cross section of the process e^e~ — > hadrons was 
the first quantity for which QCD corrections up to complete NNLO were known |90[ |9l[ . 

® DELPHI chose fif = fi and varied both within 1/2 to 2 Ecm.; ALEPH varied each independently 
within l/\/e to yTe Ecm and added both uncertainties in quadrature. Here, the uncertainties quoted by 
DELPHI are taken because they are less restrictive in limiting the range of scales. 
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For cm. energies in the e"'"e continuum, far below the TP pole and for fi = Ecm, the 
normalized cross section is given by 



a{c^e hadrons) v 2 ( 



where qi are the electrical charges of quark flavours i which are produced, i.e. for which 

Because of the above functional form of i?-y(as), relative errors of lead to absolute 
errors in Us of about the same size, AR/R ~ Aag, such that precise measurements of R 
still lead to rather large errors in ag. As another complication, electroweak corrections 
due to real Z*^ exchange and due to 7 — Z*^ interference must be applied to the data in 
the cm. energy range above 20 GeV. These corrections are negligible at Ecm — 14 GeV, 
but at Ecm =46 GeV they amount to about the same size as the QCD corrections jo^ . 

A combination of the cross section measurements of all four PETRA experiments, 
in the cm. energy range of 14 to 46.8 GeV, accounting for correlated experimental 
uncertainties, electroweak corrections and leading order initial state radiation effects, 
resulted in as(34 GeV) = 0.169 ± 0.025 [||, in NLO QCD. Applying the fuU NNLO 
QCD prediction which was not yet available at that time results in as(34 GeV) = 0.175± 
0.028 or — equivalently — in as{Mzf>) — 0.143 ± 0.018. Two groups have analysed and 
combined e^e~ hadronic cross sections measured in the cm. energy ranges from 7 to 
57 GeV Q and from 2.65 to 52 GeV Q. They both used an initially erroneous 
value of the NNLO QCD coefficient of R^, +64.8 instead of -12.8; when corrected for 
this mistake, they resuh in as(34 GeV) = 0.165 ± 0.022 ^ and in Q!s(31.6 GeV) = 
0.158 ±0.019 |||. Combining the two gives as(34 GeV) = 0.160 ±0.019 (exp. + sys.) or 
asiMzo) = 0.133 ±0.013. 

A re-analysis [H of PETRA and TRISTAN data, in the cm. energy range of 20 to 
65 GeV, took better account of higher order QED and electroweak corrections and used 
the mass of the Z" measured from LEP experiments. In NLO QCD, this study resulted 
in as{Mzo) = 0.124 ± 0.021, where the error includes experimental and systematical 
uncertainties, added in quadrature. Re-applying the NNLO QCD correction this gives 

as(42.4 GeV) = 0.175 ±0.028 or as(A/zo) = 0.126 ± 0.022 , 

which is included in the final summary section 7. A more recent determination of as 
from the total e+e~ hadronic cross section measured by CLEO at Ecm = 10.52 GeV [Q 
resulted in 

as(10.52 GeV) = 0.130 t 0:029 

and is also added to the final summary. 

On top and around the Z° resonance, the LEP experiments have collected large statis- 
tics data samples which allow accurate determinations of as- At the Z" pole, equation pT| 
must be modified according to the dominant electroweak couplings of the Z" to the 
quarks, resulting in NNLO QCD predictions of the hadronic decay width of the Z° jl^ . 
Quark mass corrections in NLO ||l^ and partly to NNLO |jl^, non-factorizable elec- 
troweak and QCD corrections |Q, top quark effects and other electroweak corrections 
apply in addition, rendering numerical expressions for Rz = T{Z'^ — > hadrons)/r(Z*' 
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leptons) rather complicated; see e.g. ref. [Q for a comprehensive report about QCD 
corrections on and Rz- 

In this review, a recent parametrisation of the NNLO QCD prediction of Rz ||l7|] , 
including all known corrections indicated above, is applied to determine as from Rz- This 
parametrisation, as given in equation is also used by the LEP collaborations in their 
combined studies of electroweak precision data . The coefficients given in equation |l5| 
are calculated for a Higgs mass Mh of 300 GeV, a top quark mass Mt of 174.1 GeV and 
for Mzo = 91.19 GeV. With the latest combined LEP result, Rz = 20.768 ± 0.024 
this gives as(Mzn) = 0.124 ± 0.004 (exp.). 

Table 2: Estimates of errors in the determination oi as{Mzo) from Rz, caused by different 
sources of uncertainties. 



error source 


Aas(Mzo) 


AMzo = ±0.0021 GeV 

AMt = ±5 GeV 

Mh = 100 ... 1000 GeV 

...4) Mza 

renormalization schemes 


±0.00003 

±0.0002 

±0.0017 

+ 0.0028 
- 0.0004 

±0.0002 


total 


+ 0.003 
- 0.002 



Errors in as{Mzo) from different sources, namely from changes of A/zo and Mt within 
their current uncertainties, from Mh in the range of 100 GeV (the current experimental 
lower limit from LEP) to 1000 GeV, and from renormalization scale uncertainties, varying 
fi from 1/4 to 4A/zo , are given in Table ^ The estimate of the scheme dependence, which 
in NNLO must be considered in addition to the renormalization scale dependence, was 
taken from reference |l99| . Neglecting the tiny errors from Mzo and from Mt, the overall 
theoretical uncertainty on as(Afzo) from Rz is estimated to be ±0.002 from Mh and Mt 
and 1 o'ooi from higher order QCD contributions. Therefore, the final result from Rz, 
to be added to the overall summary, is 

as(Mzo) = 0.124 ±0.004(exp.)±0.002(Mif,Aft)^ o.ooi(QCD) . 

Note that the precision of this result crucially depends on the assumption that the predic- 
tions of the electroweak Standard Model are strictly valid; small deviations from these 
predictions can produce large systematic shifts of as (Mzo) from Rz- The LEP data, 
however, are in excellent overall agreement with the Standard Model predictions ||l5|] , 
such that the uncertainties quoted above seem realistic. 

A combined fit to all data from LEP-I and LEP-II, including the measurement of the 
mass of the W boson, Mw, and all measured cross sections and asymmetries, instead of 
Rz alone, results in as{Mzo) = 0.120 ±0.003(exp.); a fit to all data including those from 
pp collider and lepton-nucleon scattering experiments gives as(Mzo) = 0.118±0.003(exp.) 
[ psf . These fits include simultaneous determinations of Mh and Mt, such that only the 
QCD uncertainty must be added. Because of their complicated nature, however, these 
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results are only mentioned for completeness, but are not taken to replace the above value 
of as{Mio) from Rz for the final summary. 



4.4 ds from r decays 

An important quantity to determine as from measurements at small energy scales is the 
normalized hadronic branching fraction of r lepton decays, 

r(r hadrons Vr) , . 

r(r eVeVr) 



which is predicted to be |100| 



Rr = 3.058(1.001 + (5pert + <5nonpert) ■ (23) 

Here, Jpert and ^nonpcrt are perturbative and nonperturbative QCD corrections ; (5pe rt was 
calculated to complete 0{q^) and is of similar structure to the one for Rz [BO, 100, 101|: 



i„„ . f^Hl) + 5,20 (fSa« V + 26,37 f^iiMl) V . m 



and parts of the 4* ^ ord er coefficient are also known |102]. Based on the operator product 
expansion (OPE) |l03[ , the nonperturbative correction was estimated to be small |10C|, 

iJnonpert = -0.007 ± 0.004. 

The most comprehensive determinations of as from r decays are based on recent 
studies from LEP, making use of the large data statistics available at LEP-I. The ALEPH 



[ 104 1 and the OPAL [|105| Collaborations presented measurements of the vector and the 
axial-vector contributions to the differential hadronic mass distributions of r decays, 
which allow simultaneous determination of and of the nonperturbative corrections (in 
terms of the OPE). These corrections were found to be small and to largely cancel in the 
total sum of Rr, in good agreement with the theoretical estimates. 



Table 3: Results of as{Mr) from i?r, for different variants of NNLO QCD calculations. 



Theory 


ALEPH |104| 
as{Mr) exp. theo. 


OPAL |l05| 
as{Mr) exp. theo. 


CIPT 

FOPT 

RCFT 


0.345 ±0.007 ±0.017 
0.322 ±0.005 ±0.019 


0.348 ±0.010 ±0.019 
0.324 ±0.006 ±0.013 
0.306 ±0.005 ±0.011 



The final results of asC-^r) are listed in Table ^, obtained for different variants of 
the NNLO QCD predictions: fixed order perturbation theory (FOPT) |10C||, contour 



improved perturbation theory (CIPT) |102|, expressing Jpert by contour integrals in the 
complex s-plane, and renormalon chain improved perturbation theory (RCPT) |106|, 
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where leading terms of the /?-functions are resummed by inserting so-called renormalon 
chains, i.e. gluon lines with many loop insertions. The two groups agree well on their 
as results and — approximately — on the estimated uncertainties, however different 
theoretical approaches give systematic differences in as- The FOPT results seem to 
represent the mean of these theoretical approaches. Therefore the final result from i?^, 
to be included in the final summary section 7, is taken to be 

as{Mr) = 0.323 ± 0.005(exp.) ± 0.030(theo.) , 

where the average between ALEPH and OPAL was taken and an additional error of 
±0.020, accommodating the shift between different theoretical approaches, was added in 
quadrature^ to the theoretical uncertainties given by ALEPH and OPAL. When extra- 
polated to the energy scale M^o , using the 4-loop /3- function and 3-loop matching at the 
bottom quark pole mass, this results in as{Mza) = 0.1181 ±0.0007(exp.)±0.0030(theo.). 



5 Results from hadron colliders 

Significant determinations of as from hadron collider data are obtained from bb pro- 
duction cross sections, from prompt photon production, from inclusive jet production 
cross sections and from the ratio oi W + {1 — jet) and + (0 — jet) production cross 
section, all of which arc calculated in complete NLO QCD. The latter topic, W + jet 



production, will not be discussed further here, because early measurements [108|, giving 



as{Mw) = 0.123±0.025, were put in question by new analyses showing bad disagreement 



between data and QCD, see e.g. |10!;] 



Hard scattering cross sections initiated by two hadrons with four-momenta Pi and 
P2 can be written as 

0-(A,^2)=^ jAXldx2f^{Xl,^l})fj{x2,^J.))■(T^j{pl,P2,as{^l'^),Q^/^l)) , (25) 

where pi = XiPi and p2 = X2P2 are the momenta of the interacting partons, fi and fj 
are the QCD quark and gluon distributions, Q is the characteristic scale of the hard scat- 
tering, and aij is the short-distance cross section of the hard scattering between partons 
of type i and j. This parametrization again is based — as in DIS — on the assumption 
of factorization between the short- and the long-range regimes of the scattering process, 
the transition between both being defined at the factorization scale /i / . 

In general, determinations of ag from hadron-hadron-coUisions are less precise than 
those from e^e~ annihilation or deep inelastic scattering processes, due to larger uncer- 
tainties associated with incoming hadrons: parton distributions and soft remnants from 
spectator partons, which do not participate in the hard scattering process, substantially 
add to the overall uncertainties. 



^In another study of the impact of different theoretical variants of the NNLO QCD exp ecta tion for 



Rt, it was concluded that the overall theoretical uncertainty on Q:s(Mr) is, at best, ±0.05 |107 
is larger than the respective error quoted above 



which 
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5.1 Os from bb cross sections 



The first theoretically well-defined determination of as from a purely hadronic production 



process was presented by the UAl collaboration [110|, obtained from a measurem ent o f 
the cross section of the process pp bbX for which NLO QCD predictions exist 
fe-quarks were detected through their scmilcptonic decays into muons which yield high 
transverse momenta pr with respect to the beam axis. A strong correlation between the 
decay muons and the original 6-quark allows determination of the fe-quark production 
cross section without application of a jet algorithm, thus avoiding systematic effects 
from spectator partons (or the "underlying event"). The LO QCD contribution is of 
©(ctg) and leads to a back-to-back (in azimuth) bb configuration; virtual corrections to 
this state and the emission of a third (gluon) jet are of 0{as), predicted in NLO QCD. 

Comparison of the measured cross-section for 2-body final states with NLO QCD 
predictions yielded 

as(20 GeV) - 0.1451 (exp.)+ If.l (theo.) , 

where the theoretical error includes uncertainties due to different sets of structure func- 
tions, rcnormalization/factorization scale uncertainties and the b-quark mass |11C|. 



5.2 tts from prompt photon production 

Production of high transverse momentum direct ( "prompt" ) photons in hadron collisions 
is well suited to test perturbativc QCD and to determine Qg, because photons, in contrast 
to quarks, do not hadronize and their energies and directions can — in general — be 
measured with higher accuracy than those of hadron jets. In leading order, however, 
prompt photon production is of 0{aas), compared to 0(0^) for hadron jets, and therefore 
suffers from relatively small production cross sections. In addition, there is a sizeable 
background of photons from vr*' and rj decays, such that quantitative tests of QCD from 
prompt photon production are not trivial, from an experimental point of view. 

Using complete ©(aa^) QCD predictions [112|, the UA6 collaboration determined 



from a measurement of the cross sections difference crijfp —>■ "fX) — <T{pp —> jX) |113|, 
where the poorly known contributions of the sea quarks and the gluon distributions in 
the proton cancel. The result is 

as(24.3 GeV) = 0.135 ±0.006 (exp.)l (theo.) , 

where the theoretical error includes uncertainties from the scale choice and from variation 
of the parton distribution hmctions. 



5.3 Us from inclusive jet cross sections 

The definition and reconstruction of particle jets in hadron collisions traditionally has 
followed other strategies than those used in e"'"e~ annihilation. In hadron collisions, so- 
called cone jet finders are employed which allow particles, ideally those which originate 
from the proton remnants, not to be associated with any of the reconstructed jets — 
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in contrast to the clustering algorithms used in e"'"e annihilation where all particles 
are assigned to jets; c.f. section 4.1. Nowadays, almost all of the jet studies at hadron 



collider experiments follow the "Snowmass" definition of jets |114|. Here, jets are defined 



by concentrations of transverse energy Et = \Esm9\ in cones of radius 



R = ^{Arif + , 

where t] ~ — lntan(0/2) is the pseudorapidity, (j) is the azimuthal and 9 is the polar angle 
of a particle or an energy cluster in the calorimeter of the detector, measured w.r.t. the 
point of beam crossing. 



Giele, Glover and Yu determined as |115] by fitting the NLO parton level Monte Carlo 



JETRAD |116|, which is based on the QCD matrix elements of reference [117| and the 



MRS A' parton density functions [118 , to the single inclusive jet cross sections measured 
by CDF [119 for a cone size of i? = 0.7. In each bin of Et, cts is determined for the 
scale choice = Et- The resulting values of as (St) are displayed in figure ^, where the 
error bars represent the combined statistical and theoretical uncertainties; the size of an 
additional, overall (experimental) systematic error is also indicated. The results are in 
good agreement with the QCD expectation of the running as. A corresponding overall 
QCD fit, in NLO, gives 

as(Mzo) = 0.121 ± 0.001 (stat.) ± 0.008 (sys.) ± 0.005 (theo.) ± 0.002 (pdf) 

where the theoretical error was obtained from a variation of fx between 0.5 and 2- Et, and 
the last error represents the uncertainty from using different parton density functions. 

This study was basically meant to introduce a general method and possibility to 
determine as from hadron collider jet cross sections, rather than to present a complete and 
final analysis. For instance, the parton density functions used in this study were extracted 
from (mainly deep inelastic scattering) data for a fixed input value of as(Mzo) = 0.113, 
while a coherent determination of as should allow as to vary in the density functions, 
too. Although the data statistics — and hopefully also the experimental systematic 
errors — must have improved significantly with respect to the previous data sample of 



reference [119[, and although there exist new parametrizations of density functions for 
different input values of as, no update of the measurement of as from hadron collider 
jets was published so far. Therefore the above result is retained in the final summary — 
keeping in mind, however, that the quoted uncertainties constitute lower limits rather 
than a complete assessment of the overall error. 



6 Results from heavy quarkonia decays and masses 

The mass spectra and partial decay widths of heavy quark-antiquark bound states are a 
good testing ground for QCD. For quark masses mq 3> ^qcd, the short- and long-range 
effects on the decay widths can be factorized, and the short-range part can be calculated 
by perturbative QCD. Mass splittings of heavy quarkonia states can be calculated using 
nonperturbative methods like lattice gauge theory, which also provide means to extract 
values of ag. 
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6.1 Os from quarkonia decay branching fractions 



Partial decay widths of heavy quarkonia, hke V'' = T{QQ /U+M"), TTSf and 



are calculated in NLO QCD |12C]. The expressions contain the (unknown) radial wave 
function at the origin, which can be eliminated by forming the ratios 



pggg ;L0(7r2 - 9) (^2) 



T999 ~ bajji^ 



0.4 -6.3 In 



2.6 -2.1 In 



as(/^^) 



(26) 



An early, comprehensive fit of as from the J/^* and the T branching ratios was 
presented in reference [121|, resulting in q;s(M2o) = 0.113!to;oo5. This study included 



estimates of the higher order QCD uncertainties and of relativistic corrections. 

More recently, sum rules for the T system were analysed with resummation of QCD 
Coulomb effects which are responsible for the growth of the perturbative coefficients in 
NLO QCD |l22), resulting in a^iMzo) = 0.118 ± 0.006. Similar methods were already 
applied in ealier studies [123| but have led to conflicting results. Application of NNLO 
QCD predictions resulted in |124| 



as(A'/zo) = 0.118 ±0.006 , 

which is included in the final summary of as- It should be noted that further studies of T 
sum rules in NNLO QCD |125 obtained similar results of as, however with a much more 
conservative estimate of the theoretical uncertainties — the resulting, large uncertainties 
of as may indicate that T sum rules do not appear to provide competitive results of as. 



6.2 ds from lattice calculations of mass splittings 

Early determinations of as from heavy quarkonia mass splittings were based on quenched 
approximations of lattice QCD calculations, neglecting light quark fiavour loops {Nf = 
0) . When evolved to the MS coupling at the scale M^o , these methods led to as {M^o ) = 



0.105 ± 0.004 [126|, where the error included statistical as well as estimates of systematic 
errors. 

Refined nonrelativistic QCD calculations with Nf — and Nf — 2 allowed extrapo- 
lation to the physically correct number of light quarks, Nf = 3, and included advanced 
(3-loop) perturbative extrapolation from the lattice to the MS coupling. A detailed study 
of various mass splittings of T (bb) states, which are precisely known from corresponding 



measurements, was presented in [127|. Technically, the lattice calculation for a given 
value of the "bare" lattice coupling aiat gave a value for the dimensionless quantity aS, 
where a is the lattice spacing and S is the mass splitting of suitable quarkonium states. 
From this result, divided by the measured value of E, the lattice spacing a was obtained. 
From QCD perturbation theory, the MS coupling as at the energy scale is given as 



a function of aiat 112?]. Extraction of as(a~ ) and extrapolation to as(M2o) resulted in 
as(Mzo) = 0.1174±0.0024 [ [127| , where the error included statistical as well as systematic 
uncertainties which include variations of light quark masses and estimates of truncation 
errors in the extraction of the MS coupling using perturbation theory. 
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More recently, a study based on similar lattice calculations and the same T mass 
splittings as in reference [ |127| , however using a different discretization scheme, derived 
a lower value, as{Mzo) = 0.1118 ± 0.0017 | 129[ . This is about three standard deviations 



smaller than the result from reference [127|, which led to the conclusion that the "true" 
systematic uncertainty must be three to four times larger than estimated before. 

Following this suggestion, the result on as from nonperturbative lattice calculations 
to be considered in the final summary section is chosen to be the average of the two 
results discussed above, with the overall uncertainty increased by a factor of 3, giving 

as(Mzo) =0.115 ±0.006. 



7 Summing up ... 

A summary of the as measurements discussed in the previous sections is presented in 
table H at the end of this review. The results are given, if applicable, at the relevant 
energy scale Q of the process, as{Q), and at the standard "reference" scale of the Z*^ mass, 



as(Mzo). The conversion between these two cases was done |134| using the 4-loop QCD 
expression for the running as , equation |^ with 3-loop matching at the c- and b-quark 
pole masses of 1.5 and 4.7 GeV, respectively, as discussed in section 2.5. The splitting of 
the overall uncertainties of as{Mzo), Aas(M2o), into experimental and theoretical errors 
is given in the 5*'* and the 6*'* columns of table ^ and the last column indicates the level 
of theoretical calculations on which these results are based. 

The results for as(Q), given in the 3'''' column of table ^, are presented in figure |lO| , 
together with fits of the 4-loop QCD prediction for the running as (equation |^) with 
3-loop matching at the quark pole masses. Results which were obtained from data in 
large ranges of Q are not displayed in this figure. The data are in very good agreement 
with the theoretical expectation, and prove the running of as with high significance. The 
latter point will be analysed in more detail in section 7.2. 

The values of as (A/20) are presented in figure |ll]. Within their assigned total erros, all 
results agree well with each other and with a weighted average value of as {M^o ) = 0. 1 184, 
which is indicated by the vertical line in figure 111]. The overall of all results with this 
average is 7.2 for 25 degrees of freedom, which indicates that the individual errors must 
be highly correlated. Therefore, the determination of the average value of as{Mzo) and 
its overall uncertainty requires special treatment and discussion. 



7.1 World average as(Mzo) and its overall uncertainty Aag 

The errors of most as results are dominated by theoretical uncertainties, which are es- 
timated using a variety of different methods and definitions. The significance of these 
nongaussian errors is largely unclear. Furthermore, there are large correlations between 
different results, due to common theoretical uncertainties, as e.g. for event shape mea- 
surements in e"''e"' annihilations. Correlations between as determinations from different 
processes, such as DIS and e~''e~ annihilations, or between different procedures and ob- 
servables used within the same class of processes may be present, too. Standard statistical 
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methods therefore do not apply when averaging these results. Several methods are em- 
ployed to derive an estimate of the average value 'a^{M-z,o) and its overall uncertainty, 
Aa7- The results are summarized in table 0: 



An error weighted average and an "optimized correlation" error is calculated from 
the error covariance matrix, assuming an overall correlation factor between the total 
errors of all measureme nts. This factor is adjusted so that the overall equals one 



per degree of freedom |130|. The resulting mean values, overall uncertainties and 



optimized correlation factors are given in columns 3 to 5 of Table ^, respectively. 

For illustrative purposes only, an overall error is calculated assuming that all mea- 
surements are entirely uncorrelated and all quoted errors are gaussian. The results 
are displayed in column 6. 

The simple, unweighted root mean squared of the mean values of all measurements 
is calculated and shown in column 7, labelled "simple rms" . 

Assuming that each result of as {M^o ) has a rectangular-shaped rather than a gaus- 
sian probability distribution, the resulting weights (the inverse of the square of the 
total error) are summed up in a histogram, and the resulting rms of that distribu- 



tion is quoted as "rms box" [131|. 



All of these methods have certain advantages but also include inherent problems. The 
"simple rms" indicates the scatter of all results around their common mean, but does not 
depend on the individual errors quoted for each measurement. The "box rms", which 
takes account of the errors and of their nongaussian nature, was criticized as being too 
conservative an estimate of the overall uncertainty of ag. The "optimized correlation" 
method — in the absence of a detailed knowledge of these correlations — over-simplifies 
by the assumption of one overall correlation factor. Moreover, if correlations are present, 

does not have the same mathematical and probabilistic meaning as in the case of 
uncorrelated data. In the extreme, may even be negative. 

With these reservations in mind, all four methods do provide some estimate of AoT. 
Apart from the method to calculate Aal, the result also depends on the significance 
of the data included in the averaging process: in all cases except the "uncorrelated" 
error estimate, Aa7 is largest if all data are included, and tends to smaller values if 
the averaging is restricted to results with errors Aag < Aas'""^'', i.e. if only the most 
significant results are taken into account. This is demonstrated in rows 1 to 5 of table ^. 
The dependence of Aa7 on Aai"*"^' is illustrated in figure ^2|b, where the rightmost 
results correspond to the first row of table ^ which includes all 26 as measurements 
summarized in table |[ Figure illustrates the distribution of as{Mzo) and the total 
errors of all these measurements. 

On first sight it seems logical to restrict the determination oia^{AIzo), and especially 
of Aa7, to the most significant data, if the inclusion of less significant measurements only 
enlarges AoT. Taken to the extreme, one may even be tempted to take the one result with 
the smallest quoted error as the final world average value of a7(Mzo). However, the errors 
on as(Afzo) estimated in individual studies are very often lower limits because unknown 
and additional systematic effects can only increase the total error. In some cases, small 
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Table 4: Average values of q;s(Mzo) and averaged uncertainties, for several methods to 
estimate the latter, and for several subsamples of the available data. 



row 


sample (entries) 


al(Mzo) 


opt. corr. 

Aol 


overall 

corrcl. 


uncorrel. 

Aa7 


simple rms 

Aal 


rms box 

AoT 


1 
2 
3 
4 
5 


all (26) 
Aas < 0.010 (20) 
Aas < 0.008 (18) 
Aas < 0.006 (9) 
Aas < 0.005 (4) 


0.1191 
0.1191 
0.1190 
0.1188 
0.1189 


0.0045 
0.0041 
0.0039 
0.0033 
0.0022 


0.71 
0.66 
0.62 
0.64 
0.28 


0.0012 
0.0012 
0.0012 
0.0014 
0.0017 


0.0043 
0.0037 
0.0038 
0.0029 
0.0034 


0.0057 
0.0051 
0.0050 
0.0038 
0.0033 


6 

7 
8 
9 


NNLO only (9) 
Aas < 0.008 (6) 
Aas < 0.005 (3) 
Aas < 0.004 (2) 


0.1185 
0.1184 
0.1184 
0.1175 


0.0035 

0.0031 
0.0022 
0.0026 


0.78 

0.68 
0.27 
0.95 


0.0016 

0.0016 
0.0018 
0.0019 


0.0045 

0.0026 
0.0037 
0.0006 


0.0048 

0.0032 
0.0028 
0.0019 


10 
11 
12 


only DIS (6) 
only e+e" (15) 
only pp (3) 


0.1178 
0.1209 
0.1135 


0.0040 
0.0051 
0.0074 


0.94 
0.79 
0.60 


0.0020 
0.0016 
0.0051 


0.0014 
0.0038 
0.0059 


0.0047 
0.0054 
0.0068 


13 
14 
15 


< 10 GcV (9) 
10 < ^ < 90 (9) 
g > 90 GeV (8) 


0.1177 
0.1202 
0.1213 


0.0040 
0.0064 
0.0056 


0.93 
0.56 
0.78 


0.0016 
0.0029 
0.0023 


0.0017 
0.0062 
0.0035 


0.0042 
0.0077 
0.0050 



total errors can be due to fortunate coincidences, to ignorance, over-optimism and/or 
neglect of certain error sources. In order to ensure and to test consistency of the results, 
it is therefore mandatory not to rely on a single determination alone but to possibly 
include several different, significant measurements in the averaging process. 

The averaging procedures discussed so far include results which arc based on different 
orders and types of QCD calculations, namely on lattice gauge theory and on perturbation 
theory in NLO, resummed NLO and in NNLO. Comparing and averaging these different 
types of results is nevertheless justified because all of them include estimates of the 
relevant uncertainties, such that they should be compatible with a common average 
within their given errors. In order to achieve the highest precision and confidence in a 
combined average value (^^(Afzo), however, it may be beneficial not to include results 
which are, for instance, based on calculations in lower than the maximum available order 
of perturbation theory. 
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In this sense, the averaging procedure is apphed only to those results which are based 
on complete NNLO QCD perturbation theory. Only very recently the number, the preci- 
sion and the diversity of such measurements reached a level where such a restriction still 
provides a solid basis for a meaningful averaging process which includes sufficient free- 
dom for internal consistency checks. Table ^ contains 9 measurements which are based 
on NNLO QCD; the results of averaging these are given in rows 6 to 9 of table ^, for dif- 
ferent subsamples defined by selecting measurements with a total error Aag < Aas™°^^ . 
The combined errors, Ao^, using the optimized correlation method, are displayed in 
figure [2|b. 

Requiring Aas < 0.008 rejects those as measurements which are dominated by large 
experimental uncertainties, leaving the 5 most significant NNLO results from which one 
obtains, using the optimized correlation method, 

al(Mzo) = 0.1184 ±0.0031 . 

This value is taken as the currently best estimate of the world average of as{Mio). 
According to equations ^ and ^ in 4-loop approximation and with 3-loop threshold 
matching at Mb = 4.7 GeV, this corresponds to A^^"^^ = (SlStgs) MeV and A^^^''^ = 
(29611^) MeV. 

Each of the 26 single results summarized in table ^ is compatible with this average, 
to within about one standard deviation of its assigned uncertainty or less. In order to 
investigate possible systematic deviations or trends between and within subsets of these 
measurements, the averaging procedure was repeated for deep inelastic scattering data, 
for e+e~ annihilation and for hadron collider data alone, see rows 10 to 12 of table ^, 
and for 3 energy ranges as shown in rows 13 to 15. Within the overall uncertainties (e.g. 
those derived from the optimized correlation method), all results agree well with each 
other and with the world average value derived above. Small but insignificant systematic 
differences between DIS and e"'"e~ results and between those obtained from low and from 
high energy data may be visible; these may well be accidental or may be caused, for 
example, by different methods of treating renormalization and factorization scales. 



7.2 Quantifying the running of as', determination of (3o, of Nc and 
of the functional Q-dependence 

The precision of the experimental results, the large energy range for which data are 
available, and the good agreement of the measurements with the QCD expectation of 
a running as, see figure |l^, suggest actually fitting the functional form of the energy 
dependence from the data. For this purpose, a subset of data is selected which ensures 
maximal independence between the data chosen and which is based on the most sig- 
nificant measurements in the largest possible energy range (c.f. table |^ and figure [io| ), 
namely: 

• as from r decays: as(1.778 GeV) = 0.323 ± 0.030; 

• as from moments of F2: as(2.96 GeV) = 0.252 ± 0.011; 

• as from scaling violations of F3: as(5.0 GeV) = 0.214 ± 0.021; 
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• Qfs from the hadronic width of the TP: 

• tts from hadronic event shapes: 



as(91.2 GeV) = 0.124 ± 0.005; 
as(189 GeV) = 0.110 ±0.004. 



All these results are based on complete NNLO QCD predictions — with the exception of 
the latter which includes resummed NLO calculations — and represent the most precise 
tts results of their class. The nature of their experimental and theoretical uncertainties 
ensures a maximum of independence between them. 



Table 5: Functional fits of ^^(Q^), in the energy range of 1.778 GeV < Q < 189 GeV. 
The LO QCD expectation for B = l//3o, with A^c = 3 and A/ 5, is 1.64. 



as(Q) - ... 


A 


B 


C 


^ 1 d.o.f. 


prob. 


B(i„g;)-' 




1.62 ±0.09 


0.125 ±0.032 


0.63 / 3 


0.89 


A + CIQ 


0.201 ± 0.007 
0.113 ±0.003 




(-5.3 ±0.5) • 10~4 
0.41 ±0.03 


78/3 
4.6 / 3 


0.20 



Table | summarizes some of the functional fits which were performed with these 
selected data. The general idea was to obtain qualitative measures for the functional 
form of the energy dependence of as, as well as numeric fit values for the coefficient /3o 
of the QCD /3-function, see equation H, and eventually for one of the main parameters of 
QCD, the number of colour degrees of freedom = Ca which is an integral part of /3o: 

iZTT 

Because of its simplicity and in order to avoid too many QCD-inspired biases, the func- 
tional form of the leading order QCD expression for the running as ( c.f. equation 0), 
without quark threshold matching, was fitted to the data. Higher order corrections to 
the running affect as at the smallest energy scales, Q ^ Mr, by less than 10%, which is 
well within the error of as in this energy regime. 

The results of this fit are given in the first row of table |^, with B = I/iSq = 1.62 ±0.09, 
and C = A = 0.125 ± 0.032 (in units of GeV). The of this fit is 0.63 for 3 degrees 
of freedom, which corresponds to a probability of 0.89. Within QCD, with A'c = 3 
and Nf — 5, l//3o = 1.64, which is in excellent agreement with the fit value of the 
parameter B. A = 0.125 GeV corresponds, in LO QCD, to as{Mzo) = 0.124 and to 
asiMr) — 0.31 (without threshold matching and with Nf = 5 throughout). From the 
fit value of _B = l//3o one can derive, according to equation |2^ and for Nf = 5 quark 
fiavours, 

Nc = 3.03 ±0.12 , 

which is an excellent and precise verification of the QCD group structure which implies 
Nc = 3. Alternatively, because this functional fit actually constrains the ratio B/n, 
where n is the power of Q/C inside the logarithm, the fit determines, when B is fixed to 
its QCD value of 1.64, that n = 2.03 ± 0.12. 
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Other functional forms of the energy dependence of as, which are not predicted by any 
consistent theory but which are added to demonstrate the significance of the logarithmic 
decrease of as{Q'^), are fitted and presented in the last two rows of table ||: a straight- 
line fit is clearly excluded by the unacceptable of 78 for 3 d.o.f., while a 1/Q energy 
dependence is disfavoured but cannot be excluded with the current data. 



8 Conclusions 



This review of experimental determinations of the coupling parameter of the Strong In- 
teraction, as, summarized a topical and still ongoing field of activities in the high energy 
physics community. Due to impressive theoretical developments and experimental ef- 
forts during the past 10 years, as could be determined from a large variety of physical 
observables and processes with steadily increasing precision. Current state-of-the-art 
measurements reach experimental uncertainties down to a few per cent, and theoreti- 
cal calculations which are complete up to NNLO perturbative QCD obey higher order 
uncertainties of a few per cent in Q!s(M2o), too. While not all of the available studies 
reach both these levels of precision, the most significant determinations of as, based on 
complete NNLO QCD calculations, can be summarized to a new world average value of 

a;(Mzo) = 0.1184 ± 0.0031 . 

The overall uncertainty is slightly larger than the smallest total errors quoted for some 
of the measurements, which is due to the following facts: 



• The errors of single results contain different estimates of theoretical uncertainties, 
which are not uniquely defined. 

• Many (if not all) results rely on assumptions which ultimately cannot be proven, 
like factorization between perturbative and nonperturbative effects, the nature and 
size of the latter, neglect of quark masses, exact realization of the electro-weak 
theory including the existence of a Standard Model Higgs boson, reaching of the 
continuum limit in lattice calculations etc.; the quoted theoretical uncertainties 
must therefore rather be viewed as lower limits instead of complete estimates. 

• Many of the results and their theoretical uncertainties are highly correlated, how- 
ever mostly to an unknown degree. This applies not only to their quoted errors, 
but also to the unproven assumptions mentioned under the previous item. 



A realistic determination of the world average value of as must therefore be based on 
comparing several different measurements of similar precision. It must account for possi- 
ble correlations between them and allow for underestimates and fortunate cancellations 
of errors in some cases. In the absence of any "exact" method to account for all these 
effects, the "optimized correlation" method [130 was used to calculate the world average 
ct^{Mza) and its overall remaining uncertainty as given above. The total error depends 
on the averaging method chosen - alternative procedures or preselections of data resulted 
in values ranging from ±0.0022 to ±0.0057. 
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All values of as summarized in table |^ agree well, within their quoted uncertainties, 
with the world average. Averages obtained from various subsamples of these data, like 
those from deep inelastic scattering, from e"'"e~ annihilations and from hadron colliders, 
as well as subsamples from different energy ranges, do not show significant biases or 
shifts. In fact, the measured energy dependence of as(Q^) is in excellent agreement with 
the QCD expectation and significantly proves the running of as- 

Functional fits to data in the energy range from 1.78 to 189 GeV provide evidence for 
a logarithmic (in contrast to e.g. a linear) energy dependence, and the fitted logarithmic 
slope — if interpreted in terms of Po, the LO coefficient of the QCD /?- function — provides 
an accurate value for the number of colour degrees of freedom, Nc = 3.03 ± 0.12. This is 
not actually a measurement of Nc, because many of the as determinations are based on 
theoretical predictions of observables which inherently include the QCD value of A''c = 3. 
Nevertheless, the fit result of Nc from the measured energy dependence of as constitutes 
an important consistency check between data and QCD. 

The total uncertainty of al(M2o) quoted above is 2.6%. This precision is a remarkable 
success. The error on a^{Mio), however, is much larger than those on other fundamental 
"constants" of nature, like the fine structure constant a, the weak mixing angle sin^ 9w 
or the gravitational constant Q . 

Any further reduction of the uncertainties on as will require large theoretical as well 
as experimental efforts: 



• New and higher order QCD calculations at least for some observables will be 
mandatory to understand, specify and possibly reduce theoretical uncertainties; 
most wishful candidates would be complete NNLO calculations for jet rates and 
hadronic event shapes in e"'"e~ annihilation and in DIS, which are well understood 
experimentally. 

• The precision and theoretical understanding of parton density functions and their 
correlations with extracted values of ag must be significantly improved. Future 
high statistics and high energy runs of the Tevatron hadron collider will hopefully 
provide this important input. New theoretical developments to understand and 
predict structure functions from basic principles, for instance through application 
of classical string theory, see e.g. reference |132|| , are mandatory. 

• So far, the neglect of finite (heavy) quark masses in almost all higher order QCD 
calculations potentially limits the reliability especially of results based on data 
which are close to these thresholds. While phase-space effects can be experimentally 
examined, only complete higher order QCD calculations for massive quarks, see e.g. 
reference |9), can close this gap of potentially large systematic uncertainties. 

• The role of nonperturbative, long distance effects and the application of factoriza- 
tion in many processes and as determinations must be further investigated and 
understood. Power corrections replacing or supplementing hadronization models, 
see e.g. references k3, b5|, or higher twist corrections from lattice QCD, see e.g. 



reference |133], are promising new developments which have the potential for new, 



significant insights and improvements. 
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The physics of hadronic interactions at high energy colhciers is rich and colourful. 
Much has been achieved in the past, however there are still many fundamental and open 
questions. Further and more precise determinations of the strong coupling parameter 
will continue to concern and employ many motivated scientists in the future. 
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Figure 2: (a) The running of as{Q), aceording to equation in 1-, 2- and 3-loop ap- 
proximation, for Nf = 5 and with an identical value of = 0.22 GeV. The 4-loop 
prediction is indistinguishable from the 3-loop curve, (b) Fractional difference between 
the 4-loop and the 1-, 2- and 3-loop presentations of as{Q), for Nf = 5 and chosen 
such that, in each order, as(Afzo) = 0.119. 
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Figure 3: (a) 4-loop running of as{Q) with 3-loop quark threshold matching according 
to equations 0and||, with A.|^~^'' = 220 MeV and charm- and bottom-quark thresholds 

at the pole masses, fj.i'^^^^'^ = = 1.5 GeV and = Ah ^ 4.7 GeV (fuU line), 

compared with the unmatched 4-loop result (dashed line), (b) The fractional difference 
between the two curves in (a). 
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Figure 4: (a) asiM^o) determined from the scaled hadronic width of the Z°, i?z — 20.768, 
in leading, next-to-leading and in next-to-next-to leading order QCD, as a function of the 
renormalization scale factor — fi/Mio- NLO solutions according to the PMS (•), the 
EC (o) and the BLM (o) scale optimization methods are marked, (b) Scale dependence 
of the QCD coefficients i?„ of Rz- 
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Figure 5: Comparison of CCFR data M with NNLO QCD fits P§ 
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Figure 6: xF^ as a function oi x at — 5 GeV^, with a power law fit to the region 
X < 0.1 (left figure; logarithmic x-scale) and the shape of the SLAG F2 data for x > 0.5 
(right figure; linear x-scale) which are used to calculate the GLS integral in the shaded 
regions (figures from p6|). 
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Figure 8: Running of as from hadronic event shapes at LEP, measured by L3. The 
results at energies below 91 GeV are from radiative events at 2Ebeam ~ M^o (figure from 
reference 
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Figure 9: Values of as{ET) extracted from CDF inclusive jet cross section data (figure 
adapted from [115|). 
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Figure 10: Summary of as{Q) 
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Figure 11: Summary of as(Af2o) (fiilcd symbols represent results based on complete 
NNLO QCD). 
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Figure 12: (a) Results of as{Mzo) versus their respective overall errors; full symbols are 
based on complete NNLO QCD. (b) Dependence of A57 on the selection of results with 
errors Aasl^^z") ^ Aai"^°'^\ for different methods to calculate a7(Mzo). The rightmost 
data points correspond to no cut on Aas™"^'* . 
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Table 6: World summary of measurements of ag (DIS = deep inelastic scattering; 
GLS-SR = Gross-Llewellyn-Smith sum rule; Bj-SR = Bjorken sum rule; (N)NLO = 
(next-to-)ncxt-to-leading order perturbation theory; LOT = lattice gauge theory; resum 
= resummed NLO). 
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